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Catechol-O-methyltransferase (COMT) is an enzyme that is particularly important for the 
metabolism of dopamine. Functional polymorphisms of COMT have been implicated 
in working memory and numerical cognition. This is an exploratory study that aims 
at investigating associations between COMT polymorphisms, working memory, and 
numerical cognition. Elementary school children from 2th to 6th grades were divided into 
two groups according to their COMT val158met polymorphism [homozygous for valine 
allele (n — 61 ) vs. heterozygous plus methionine homozygous children or met+ group (n = 
94)]. Both groups were matched for age and intelligence. Working memory was assessed 
through digit span and Corsi blocks. Symbolic numerical processing was assessed through 
transcoding and single-digit word problem tasks. Non-symbolic magnitude comparison 
and estimation tasks were used to assess number sense. Between-group differences 
were found in symbolic and non-symbolic numerical tasks, but not in working memory 
tasks. Children in the met+ group showed better performance in all numerical tasks 
while val homozygous children presented slower development of non-symbolic magnitude 
representations. These results suggest COMT-related dopaminergic modulation may be 
related not only to working memory, as found in previous studies, but also to the 
development of magnitude processing and magnitude representations. 

Keywords: catechol-O-methyltransferase, COMT, numerical cognition, non-symbolic magnitude, symbolic 
magnitude, dopamine 



INTRODUCTION 

Dopamine, a neurotransmitter of the monoaminergic cate- 
cholamine group, is engaged in several brain and cognitive 
functions. Traditionally, dopamine is implicated in informa- 
tion processing in prefrontal-related working memory tasks 
(Dickinson and Elvevag, 2009; Mier et al., 2010). Among other 
domains, working memory is important for numerical cogni- 
tion (Geary, 1993; Rotzer et al, 2009; Geary et al, 2012). For 
this reason, dopamine determines performance in numerical 
tasks always when working memory is required. Besides gen- 
eral cognitive factors such as working memory, specific cogni- 
tive capacities such as symbolic and non-symbolic magnitude 
representation also determine mathematics abilities (Mazzocco 
et al., 2011). Interestingly, there is compelling evidence relat- 
ing basic counting processes to dopamine concentration on the 
prefrontal cortex (PFC; Allman et al, 2012). Could it then 
be that a more direct link between dopamine function and 
math achievement also exists? In case there is a direct con- 
nection between dopamine and arithmetic performance, which 
aspects of numerical processing are more influenced by vari- 
ability in dopamine bioavailability? In this article, we set out 
to explore Catechol-O-methyltransferase (COMT) variability 



in the population as a clue to disentangle a hypothetical 
connection between dopamine bioavailability and numerical 
cognition. 

COMT is important to the degradation of several cathe- 
cholamines, but it is especially relevant to the metabolism of 
dopamine, as it is responsible for degrading more than 60% of 
PFC dopamine (Karoum et al., 1994). Thus, COMT polymor- 
phisms allow investigation of interactions between dopaminergic 
and cognitive functions. The COMT is one of the enzymes 
that inactivate dopamine in the synaptic cleft (Standaert and 
Galanter, 2009). A nucleotide substitution (G to A) at codon 
158 of the COMT gene creates a functional polymorphic vari- 
ation (valinel58methionine, or just vall58met) of this enzyme, 
which is especially important in the PFC (Weinshilboum et al., 
1999). According to Chen et al. (2004), the enzyme metl58met 
has a degradation rate between 1/3 and V4slower than the enzyme 
vall58val. Individuals who are homozygous for the valine allele 
(val) have higher COMT activity and lower concentrations of 
extracellular dopamine. Furthermore, individuals homozygous 
for the methionine allele (met) tend to have lower levels of COMT 
activity and, consequently, higher bioavailability of dopamine 
in the synaptic cleft. Also according to Chen et al. (2004), the 



www.frontiersin.org 



August 2013 | Volume 4 | Article 531 | 1 



Julio-Costa et al. 



COMT and numerical cognition 



heterozygous polymorphism (vall58met) has an intermediate 
enzyme activity. 

Different patterns of information processing have been asso- 
ciated with the vall58met COMT polymorphism (Stein et al., 
2006). These differences are assumed to reflect variability in 
dopamine bioavailability at the synaptic cleft in the PFC. 
Homozygous metl58met individuals tend to have what is called 
in the literature a more "cautious style" of information process- 
ing, with better focusing of attention as well as higher work- 
ing memory capacity ("worrier strategy," Stein et al., 2006). 
Otherwise, vall58val homozygous individuals are prone to having 
a more exploratory style, which favors impulsivity and adapta- 
tion under stress conditions ("warrior strategy"). A meta-analysis 
of adult studies showed that there is an association between the 
COMT genotype and PFC activation, with a moderate to large 
effect size (d = 0.73, without evidence for publication bias; Mier 
et al., 2010). This study also showed that, in addition to exhibiting 
a lower PFC activation, the methionine allele carriers performed 
better in working memory tasks, while the valine allele carriers 
exhibited higher PFC activation levels and worse behavioral per- 
formance in these tasks. Although some studies fail to disclose 
differences between the COMT genotypes on working memory 
tasks (MEls et al, 2004; Tan et al, 2007; Blanchard et al, 2011; 
Wardle et al., 2013), meta-analytic evidence is compelling and less 
biased than the results of single studies. 

Research based on child samples is scarce. Some studies sup- 
port the view that met carriers are favored in working memory 
tasks, such as demonstrated by Diamond et al. (2004) in chil- 
dren aged between 6 and 14 years old. In addition, COMT has 
been found to modulate adolescents' behavior on two tasks of 
working memory (digit span and Corsi blocks), but not on a 
motor speed measure (Wahlstrom et al., 2007). Other evidence 
suggests that the association between the COMT polymorphisms 
and working memory may be moderated by age. Thus, analyz- 
ing individuals aged between 6 and 20 years old, Dumontheil 
et al. (2011) found that polymorphic-related working memory 
differences were present only in children older than 10 years. 

The relevance of working memory for mathematical learning 
is well-established, regarding both typical and atypical develop- 
ment (Raghubar et al., 2010). As an example, in previous research 
(Costa et al., 2011), we found that working memory effects on 
simple word problems were attenuated by including finger gnosia 
as covariate in the analyses. As finger discrimination abilities pre- 
sumably enable finger counting, this result may be interpreted 
in terms of the offloading mechanisms proposed by Alibali and 
DiRusso (1999, see also Costa et al., 2011). The executive com- 
ponent of working memory is more strongly implicated in math 
performance than storing mechanisms (Raghubar et al., 2010). 

Another important contribution of working memory is related 
to number transcoding (changing the type of numerical represen- 
tation, for example: "4" — > "four"). Children with higher count- 
ing span show superior performance in a numerical transcoding 
task compared to those with lower span (Camos, 2008). In 
that study, working memory level was also correlated with the 
complexity of the items to be transcoded. Another study exam- 
ining the transcoding abilities of 7-years-old children suggested 
that visuospatial working memory predicts the performance in 



transcoding tasks in the context of German language, in which 
there is decade-unit inversion in two-digit numbers (Zuber et al., 
2009). According to Barrouillet et al. (2004), working memory 
resources are increasingly required for transcoding complex num- 
bers. Evidence also indicates that transcoding may be performed 
as a purely algorithmic procedure, not requiring access to seman- 
tic representations of magnitude (Deloche and Seron, 1982a,b; 
Barrouillet et al, 2004). 

The connection between the vall58met COMT polymor- 
phism, working memory, and mathematical achievement was 
explored in a study of MRI-related events by Tan et al. (2007). 
Individuals who carried the valine allele showed higher activa- 
tion levels of the dorsolateral PFC in comparison to individuals 
with the met/met COMT genotype even in the absence of dif- 
ferences in behavioral performance. COMT activity levels were 
positively correlated with the maintenance and processing of 
information in working memory during mathematical tasks, but 
not with retrieval operations from long-term memory. These 
pieces of evidence indicate that individuals with lower dopamine 
levels (met/met homozygous) can perform the same operations 
presumably recruiting less working memory resources. 

Besides working memory capacity, several other abilities and 
cognitive representations such as magnitude representation are 
recruited in number processing and calculation (Prado et al., 
2011). The non-symbolic magnitude representation is the capac- 
ity to analogically, intuitively, and automatically discriminate, 
estimate, and compare magnitudes. The non-symbolic magni- 
tude representation can be observed already in infants (Xu et al., 
2005) and is very basic (Dehaene, 1992). Operations with non- 
symbolic magnitudes are executed in an approximate manner. 
Accuracy of the magnitude representation may be assessed by the 
Weber fraction (the minimal proportional numerical difference 
that may be discriminated between two sets of objects), or by the 
coefficient of variation in estimation of the numerosity of sets 
(Gilmore et al., 2011). This inverse association between numer- 
ical distance and reaction times and error rates is known as the 
distance effect, and it represents an instantiation of the Weber law 
(Dehaene, 1989; Dehaene et al., 1990). 

A direct link between dopaminergic modulation and magni- 
tude estimation is postulated in the accumulator model (Meek 
and Church, 1983). This model comes from temporal informa- 
tion processing research (Gibbon, 1977), and was successfully 
generalized to numerical cognition (Meek and Church, 1983; 
Dormal and Pesenti, 2012). It proposes the existence of a pace- 
maker that emits pulses which are gated into an accumulator. The 
pulse transmission is regulated by two different modes: a counter, 
when operating in the event mode (numerical estimation) or a 
clock, when operating in the run mode (time estimation). The 
accumulator value can be held in working memory and compared 
to a previous accumulator value stored in reference memory. A 
decision process determines the appropriate response by discrim- 
inating which of the two values is the larger one (Meek and 
Church, 1983; Allman et al, 2012). Meek and Church (1983) 
obtained evidence in rodents that learning of numerosity and 
time estimation have similar psychophysical functions, interfere 
one with another and are similarly affected by pharmacological 
manipulations. These results suggest that a common magnitude 
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processing mechanism underlies time and numerosity estima- 
tions (Walsh, 2003; Dormal and Pesenti, 2012). 

Meek (2006) showed, in animals that lesions in frontal areas 
led to selective reduction of the modulatory control of the clock 
speed, interfering with estimation of duration. Results also indi- 
cated that this speed reduction was mediated by dopamine recep- 
tors. More precisely, dopamine agonists have a dose-dependent 
influence on magnitude estimation (Meek and Church, 1987; 
Lustig and Meek, 2005; Cheng et al., 2008). It was veri- 
fied that the administration of methamphetamine (a dopamine 
agonist) produced an overestimation in time and numerical 
responses of rats. In other words, higher dopamine bioavail- 
ability (provocated by the methamphetamine) was associated 
with a distorted perception of time and numbers in a non- 
human sample (Meek and Church, 1983). Balci et al. (2013) 
showed that the COMT polymorphisms that lead to higher pre- 
frontal dopamine resulted in weaker manifestation of memory 
variability in timed behavior. More specifically, the associa- 
tion between COMT and time comparison has been studied 
by Wiener et al. (2011), who showed that the COMT poly- 
morphism influences the perception of time magnitude, with 
individuals with at least one methionine allele presenting a 
better performance (smaller coefficients of variation). Thus, 
it is plausible to establish a relationship between COMT and 
the mechanisms of magnitude representation (Coull et al., 
2011). 

The role of working memory and magnitude representa- 
tions in arithmetic learning suggests, therefore that an important 
source of interindividual variation may be related to COMT 
polymorphisms. They can have an impact on dopaminergic func- 
tions related to working memory and also possibly to magnitude 
representations. This exploratory study aims at: (1) investigat- 
ing the influence of dopamine bioavailability, operationalized by 
means of vall58met polymorphism of the COMT enzyme, on the 
performance of working memory tasks; (2) assessing possible dif- 
ferences between COMT polymorphic groups on numerical tasks; 
and (3) exploring the interplay between dopamine bioavailabil- 
ity, working memory, and symbolic and non-symbolic numerical 
tasks. We hypothesize that met carriers will have better perfor- 
mance on numerical tasks requiring working memory abilities 
due to their higher levels of dopamine in the synaptic cleft. 
The relevance of COMT polymorphisms for non-symbolic rep- 
resentations is still speculative. However, based on the previously 
presented evidence, one can argue that met carriers should also 
have a better performance on basic non-symbolic numerical 
tasks. Therefore, this study also aims at investigating whether 
there are group differences in non-symbolic tasks between the 
different COMT allele combinations. If the effects of COMT 
on number processing are not completely mediated by work- 
ing memory capacity, between group differences should remain 
present in numerical tasks even after removing the impact of 
working memory capacity. 

MATERIALS AND METHODS 
PARTICIPANTS 

Two hundred and nineteen school-aged children were recruited 
from public and private schools from Belo Horizonte, Brazil. 



This study was approved by the local research ethics commit- 
tee. Participation occurred only after informed consent was 
obtained in written form from parents or the legal representatives, 
and orally from children. All children had normal intelligence 
and did not present any major psychiatric or neurological ill- 
nesses. Data from 64 children were not included in the analyses, 
because they did not complete one of the tests (n = 33), had 
a poor R 2 on the fitting procedure to calculate their internal 
Weber fraction on the non-symbolic comparison task (R 2 < 
0.2; n = 13) and/or were outliers and extreme cases on any 
of the instruments used (n = 18), which might suggest they 
were not able to solve the task or were not paying adequate 
attention to it. 

The final sample was constituted by 155 children with ages 
ranging from 8 to 12 years. They were attending elementary 
school grades ranging from 1st to 6th. Sixty-one children were 
homozygous for the valine allele (val/val), 18 children were 
homozygous for the methionine allele (met/met), and 76 chil- 
dren were heterozygous (val/met). No deviations from Hardy- 
Weinberg Equilibrium were observed in our study (x 2 = 0.001; 
p > 0.05). Due to the small number of homozygous subjects for 
the methionine allele, we divided the sample into two broader 
groups: (a) children homozygous for the valine allele (thereafter 
val/val, n = 61); and, (b) children with at least one methionine 
allele (thereafter met+, n = 94). These two groups were matched 
by age, grade, sex and intelligence (allp > 0.05). Only 7% of the 
final sample was left-handed. Demographic data are described in 
Table 1. 

GENETIC ANALYSIS 

Ten milliliters of peripheral blood samples were collected and 
DNA was extracted with proteinase K digestion followed by a 
saline precipitation protocol (Miller et al., 1988). DNA samples 
were quantified using an optical absorbance spectrophotome- 
ter. COMT vall58met polymorphism was genotyped with tetra- 
primer amplification refractory mutation system-polymerase 
chain reaction (ARMS-PCR), as described in the literature (Ruiz- 
Sanz et al, 2007). ARMS-PCR reactions contained 200-400 ng 
of genomic DNA, 0.1 u,M of PI primer, 0.05 [iM of P2 primer, 
0.75 [iM of P3 primer, 0.1 [iM of P4 primer, 0.5 (iM of dNTPs, 
5% DMSO, and 0.1 U Taq DNA polymerase in a 20 |xL final vol- 
ume. The PCR cycling program included an initial denaturation 
step at 94°C for 4min, followed by 30 cycles of denaturation 
at 94°C for 30s, annealing at 58°C for 30s, and extension at 



Table 1 | Descriptive data of the sample distribution according to 
COMT genotype. 




val/val 


met+ 




df 


P 


N 

Sex (%female) 


61 
59.00 


94 
59.60 


<0.01 


1 


0.95 




Mean (SD) 


Mean (SD) 


t 


df 


P 


Age (years) 
Raven (z-score) 


9.61 (1.20) 
0.43 (0.75) 


9.81 (1.10) 
0.52 (0.70) 


-1.08 
-0.82 


153 
153 


0.28 
0.41 
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72° C for 45 s, which was followed by last extension step extension 
at 72°C for 5 min. ARMS-PCR products were resolved elec- 
trophoretically on 6% polyacrylamide gels, stained with silver 
nitrate. This ARMS-PCR amplification system may produce three 
amplicons. A 626 bp amplicon is produced by amplification of P 1 
and P2 primer pair, functions as an intratube PCR control, and 
is expected to amplify in every sample containing human DNA. 
A 222 bp amplicon is produced by the amplification with PI and 
P4 primer pair and corresponds the met allele. A 45 1 bp ampli- 
con is produced by the amplification with P2 and P3 primer pair 
corresponds to the val allele. Approximately, 20% of the ampli- 
fication reactions were confirmed by restriction with the enzyme 
Hsp92ll. 

PROCEDURE 

The data collection occurred in the children's schools and was 
divided in three distinct phases. The first phase was conducted 
in group. Children were tested using an intelligence test (Raven's 
Colored Progressive Matrices) and the Arabic number writing 
task. Afterwards, parents or the legal representatives were invited 
to their children's schools. In this occasion, biological sample 
collection (peripheral blood) was performed. The third phase 
was the individual neuropsychological assessment. Children were 
tested in a quiet room in their own schools using paper and 
pencil neuropsychological tasks, the only exceptions being the 
non-symbolic numerical tasks, which were presented in a 15.4 
inches laptop. 

INSTRUMENTS OF NEUROPSYCHOLOGICAL ASSESSMENT 

Raven's colored progressive matrices 

General intelligence was assessed with the age-appropriate 
Brazilian version of Raven's Colored Matrices (Angelini et al., 

1999) . The z-scores were calculated based on the manual's 
norms. 

Digit span 

Verbal short-term and working memory were assessed with 
the Brazilian WISC-III Digits subtest (Figueiredo, 2002). 
Performance in the forward order was considered a measure 
of phonological short-term memory and the backward order 
was used to assess verbal working memory (Figueiredo and 
Nascimento, 2007). We evaluated the total score (correct trials x 
span) in the forward and the backward orders (Kessels et al., 

2000) . 

Corsi blocks 

This test is a measure of the visuospatial component of 
short-term and working memory. It is constituted by a set 
of nine blocks which are tapped, in a certain sequence, by 
the examiner and this child must repeat the same sequence. 
The test starts with sequences of two blocks and can reach 
a maximum of nine blocks. We used the forward and back- 
ward Corsi span tasks according to Kessels et al. (2000). 
In the forward condition, the child is instructed to tap the 
blocks on the same order as the examiner, in the backward 
condition, in the inverse order. We evaluated span in each 
condition. 



Arithmetic word problems 

Twelve arithmetical word problems were presented to the child on 
a sheet of paper while the examiner read them aloud simultane- 
ously to avoid reading proficiency bias. There were six addition 
and six subtraction items, all of them with single-digit operands 
and results ranging from 2 to 9 (e.g., "Annelise has 9 cents. She 
gives 3 to Andressa. How many cents does Annelise have now?"). 
The child had to solve the problems mentally and write the answer 
down in the Arabic format as quickly as possible. Cronbach's a 
of this task was 0.83 (Costa et al., 2011). We computed the total 
correct response in the task. 

Arabic number writing task (Verbal to Arabic transcoding task) 

This is a symbolic task on which children were instructed to write 
the Arabic form of dictated numbers. This task is constituted by 
40 items, up to 4 digits (3 one-digit numbers, 9 two-digit num- 
bers, 10 three-digit numbers, and 18 four-digit numbers). The 
percentage of the correct responses in the task was computed. 
The internal consistency of this task is 0.96 (KR-20) (Lopes-Silva, 
2013). 

Arabic number reading task (Arabic to verbal transcoding task) 

Twenty eight Arabic numbers printed in a booklet were presented 
one at a time, to the children, who were instructed to read them 
aloud. The item set consists of numbers up to 4 digits (3 one- 
digit numbers, 9 two-digit numbers, 8 three-digit numbers, and 8 
four-digit numbers). The total correct response was computed. 
The internal consistency of the task was 0.90 (KR-20 formula) 
(Lopes-Silva, 2013). 

Non-symbolic magnitude estimation task 

In the non-symbolic magnitude estimation task, participants 
were asked to estimate with a verbal response the quantity of dots 
presented on the computer screen. Black dots were presented on 
a white circle over a black background. Numerosities were 1, 2, 3, 
4, 5, 10, 16, 24, 32, 48, 56, or 64 dots. Each numerosity was pre- 
sented 5 times, each time in a different configuration and in such 
a way that the same numerosity never appeared in consecutive tri- 
als. The task comprised 60 testing trials. The maximum stimulus 
presentation time was 1000 ms, fast enough to avoid counting. 
Inter-trial interval was 700 ms. As soon as the child responded, 
the examiner, who was sitting next to the child, pressed the space- 
bar on the keyboard and typed the child's answer. Between each 
trial, a fixation point appeared on the screen for 500 ms — a cross, 
printed in white, with 3 cm in each line (see Figure 1 A). As a 
measure of non-symbolic magnitude representation acuity, we 
calculated the mean coefficient of variation (estimation cv mean) 
of the numbers ranging from 10 to 64 of the responses for each 
child (Pinheiro-Chagas, 2012). Numbers between 1 and 5 were 
not included in the analyses because they are on the subitizing 
range, which require an exact access to non-symbolic magnitude 
and it is not on the scope of the present work. 

Non-symbolic magnitude comparison task 

In the non-symbolic magnitude comparison task, the participants 
were instructed to compare two simultaneously presented sets of 
dots, indicating which one contained the larger number. Black 
dots were presented on a white circle over a black background. 
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On each trial, one of the two white circles contained 32 dots (ref- 
erence numerosity) and the other one contained 20, 23, 26, 29, 35, 
38, 41, or 44 dots. Each magnitude of dot sets was presented eight 
times. The task comprised 8 learning trials and 64 experimental 
trials. Perceptual variables were varied such that in half of the tri- 
als individual dot size was held constant, while in the other half 
the size of the area occupied by the dots was held constant (see 
exact procedure descriptions in Dehaene et al., 2005). Maximum 
stimulus presentation time was 4000 ms, and inter-trial interval 
was 700 ms. before each trial, a fixation point appeared on the 
screen — a cross, printed in white, with 30 mm in each line. If the 
child judged that the right circle presented more dots, a prede- 
fined key localized in the right side of the keyboard should be 
pressed with the right hand. On the contrary, if the child judged 
that the left circle contained more dots, then, a predefined key 
on the left side had to be pressed with the left hand (Costa et al., 
2011; see Figure IB). As a measure of the acuity of the ANS, the 
internal Weber fraction (thereafter w) was calculated for each child 
individually. The proportion of responses in which the target was 
considered to be larger than the reference was determined for 
each numerical ratio tested. The proportion of these responses 
obtained for each ration were plotted and a Log-Gaussian pre- 
dictor was estimated. The R 2 obtained for each child was used 
as a reference for the individual goodness-of-fit. Children with 
R 2 < 0.2 were excluded from all further analyses. The slope of the 
Log-Gaussian corresponds to w (Dehaene, 2007), and describes 
the amount by which one would need to change a particular target 




FIGURE 1 | (A) Non-symbolic magnitude estimation task. (B) Non-symbolic 
magnitude comparison task. 



in order for a given participant to correctly detect that it is larger 
than the reference on 75% of the trials (Piazza and Dehaene, 
2004). 

RESULTS 

Analyses of variance (ANOVA) were run to identify significant 
group differences between met+ and val/val groups regard- 
ing working memory and numerical abilities. To investigate the 
extent to which the Weber fraction (w) and the mean coef- 
ficient of variation from the non-symbolic magnitude estima- 
tion task (cv mean) contributed to these group differences, 
ANCOVA models including these two variables as covariates 
were calculated. Analyses in which these non-symbolic mag- 
nitude processing tasks removed group differences were inter- 
preted as indicative of a specific influence of non-symbolic 
magnitude representations on other numerical and arithmetic 
processes. 

Even though there was no statistical difference between groups 
regarding age (p = 0.28), after calculating the effect size of 
this difference (<i = 0.18) we included age as a covariate in 
all further analyses. Generally, all results remained qualitatively 
unchanged, independently of age being statistically controlled 
or not. Therefore, we decided to focus only on the one-factor 
ANOVA. We also examined the correlations between all numer- 
ical tasks for each of the allelic groups separately. We considered 
p-values below 0.05 as being significant. 

COMT POLYMORPHISMS AND WORKING MEMORY 

No group difference between val/val and met+ groups was 
observed in working memory tasks (p > 0.15). We compared the 
groups with a series of ANOVAs on Corsi Blocks and Digits Span 
and all F tests were lower than 1.00 [with the exception of Digit 
Span forward, on which -F(i53) = 2.12] and all the effect sizes were 
significantly small {x\ 2 < 0.01). 

COMT POLYMORPHISMS AND NUMERICAL AND ARITHMETIC TASKS 

Our second main goal was to investigate if there were any group 
differences in numerical tasks. Results of group comparisons are 
depicted in Table 2. 

Statistically significant differences between val/val and met+ 
were found in all numerical tasks. Children in the met+ group 
presented superior performance in all numerical and arithmetical 
tasks (Table 2). After controlling for the impact of w in ANCOVA 
models, all comparisons remained significant, with the exception 



Table 2 | Analysis of variance (ANOVA) and covariance (ANCOVA) of the numerical tasks. 



Tasks 


val/val (n = 61) 


met+ (n = 94) 




ANOVA 




ANCOVA 


(Covariate: Weber fraction) 


Mean (SD) 


Mean (SD) 


F 


P 


I 2 


F 


P 




Arithmetic word problems 


9.15 (2.57) 


9.91 (1.98) 


4.38 


0.04 


0.03 


2.97 


0.09 


0.02 


Arabic number reading task 


25.97 (3.15) 


27.15 (1.75) 


8.98 


<0.01 


0.06 


6.94 


0.01 


0.04 


Arabic number writing task 


0.85 (0.22) 


0.95 (0.08) 


15.64 


<0.01 


0.09 


12.11 


<0.01 


0.07 


Non-symbolic magnitude 


0.31 (0.16) 


0.26 (0.10) 


6.13 


0.01 


0.04 








comparison task (w) 


















Non-symbolic magnitude 


0.19(0.07) 


0.17 (0.06) 


4.65 


0.03 


0.03 


1.96 


0.16 


0.01 



estimation task (cv mean) 
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of arithmetic word problems and estimation cv mean (Table 2). 
Although not significant, it is important to note that the effect 
sizes associated with these two comparisons were moderate or 
small and that after removing the effect of the non-symbolic mag- 
nitude estimation (cv mean) the effect sizes remained roughly 
unchanged. 

Moreover, the between-group differences in arithmetic word 
problems and w disappeared [F(is2) = 2.50, p = 0.12, r| 2 = 0.02; 
F(i52) = 3.41, p = 0.07, r| 2 = 0.02, respectively] when cv mean 
was the covariate. One can infer that both w and cv mean are 
functionally associated: the difference on w disappears when cv 
mean is used as a covariate and vice versa. 

Finally, group differences in the Arabic number reading and 
Arabic number writing tasks also remained significant [F(i52) = 
6.13, p = 0.01, r| 2 = 0.04; F (152) = 11.45, p < 0.001, r| 2 = 0.07, 
respectively] when cv mean was the covariate. This suggests 
that differences in transcoding skills between val/val and met+ 
cannot be attributed to the effect of non-symbolic magnitude 
representations. 

To further explore the association between COMT poly- 
morphisms, working memory and numerical achievement, we 
investigated the correlation between these variables in each 
group separately. The section above the diagonal refers to the 
met+ group (gray background) whereas the one below the 
diagonal concerns the correlation values of the val/val group 
(Table 3). Significant correlations were found in both groups. 
Arabic number reading, Arabic number writing and arithmetic 
word problems correlated with each other (Table 3). Moreover, 
digit span backward correlated with Arabic number reading 
task, Arabic number writing task and non-symbolic magni- 
tude estimation task. Corsi blocks backward presented signifi- 
cant correlations to Arabic number writing task. Regarding the 
non-symbolic tasks, the magnitude estimation task correlated 
with the magnitude comparison task and with arithmetic word 
problems. 

Similar correlation patterns were observed in the val/val and 
met+ groups. At first, the correlation patterns found in met 



carriers were analyzed. Inspection of Table 3 (upper diagonal) 
reveals that digit span backward did not correlate with the non- 
symbolic magnitude comparison task. All of numerical tasks 
correlated with each other, with the exception of Arabic number 
reading that did not correlate to w nor to estimation cv mean and 
the Arabic number writing task did not correlate to estimation 
cv mean. Surprisingly, there was a significant correlation between 
Arabic number writing task and w. 

Children with the val/val genotype presented a different pat- 
tern of correlations (lower diagonal of Table 3). Arabic number 
reading, Arabic number writing and arithmetic word prob- 
lems tasks presented significant correlations with non-symbolic 
magnitude estimation task, but did not correlate with the non- 
symbolic magnitude comparison task. All of these symbolic 
numerical tasks also correlated with digit span and Corsi blocks 
(backward orders), with except for the arithmetic word problems. 

In summary, in both groups, Arabic number reading, Arabic 
number writing, and arithmetic word problems that rely more 
strongly on working memory, presented significant correlations 
with digit span and Corsi blocks backward. In general, both 
val/val and met+ groups presented a very distinctive pattern of 
correlations and only in the met+ group all of these numerical 
tasks were correlated to the non-symbolic ones. 

DISCUSSION 

In the present study the association between COMT genotypes 
in vall58met polymorphism and numerical cognition was exam- 
ined in school-aged children. Our results can be summarize as 
follows: (1) Children with at least one met allele (met+) showed 
better performance in numerical tasks when compared to chil- 
dren with val/val genotype. Numerical tasks encompassed simple 
operations (arithmetical word problems), transcoding (Arabic 
number reading and Arabic number writing tasks) and non- 
symbolic magnitude processing. Effect sizes of these difference 
were moderate or small, but these results were in accordance to 
what was expected, since the groups were stratified using a single 
gene as criterion; (2) when Weber fraction (w) was entered as a 



Table 3 | Correlations between the neuropsychological measures of met group in the upper diagonal (gray back ground) and val/val group in 
the lower diagonal (white background). 



123456789 10 



1. Raven (z-score) 


1 


-0.03 


0.23* 


0.07 


0.30** 


-0.15 


0.02 


0.12 


-0.04 


-0.09 


2. Digit span (forward) 


0.15 


1 


0.27* 


0.06 


0.23* 


0.07 


0.27** 


0.39** 


-0.23* 


-0.02 


3. Digit span (backward) 


0.36»* 


0.38** 


1 


-0.06 


0.29** 


0.22* 


0.31** 


0.42** 


-0.13 


-0.24* 


4. Corsi blocks (forward) 


0.17 


0.21 


0.22 


1 


0.22* 


0.70 


0.07 


0.07 


0.04 


0.07 


5. Corsi blocks (backward) 


-0.01 


0.04 


0.20 


0.41** 


1 


0.76 


0.25* 


0.33** 


-0.24* 


-0.26* 


6. Arabic number reading task 


0.27* 


0.15 


0.29* 


0.31* 


0.25* 


1 


0.61** 


0.34** 


-0.12 


-0.05 


7. Arabic number writing task 


0.17 


0.26* 


0.34** 


0.24 


0.38** 


0.79** 


1 


0.49** 


-0.29** 


-0.09 


8. Arithmetic word problems 


0.19 


0.21 


0.23 


0.09 


0.14 


0.49** 


0.53** 


1 


-0.31 


-0.24** 


9. Non-symbolic magnitude 


-0.09 


-0.23 


-0.24 


-0.23 


-0. 14 


-0.14 


-0.76 


-0.07 


1 


0.35** 


comparison task (w) 

10. Non-symbolic magnitude 

estimation task (cv mean) 


-0.28 


-0.21 


-0.37** 


-0.11 


-0.06 


-0.42** 


-0.46** 


-0.25* 


0.31* 


1 



*p < 0.05; **p < 0.07. Correlations shown in italics were significant in the group being analyzed, but not in the other one, while correlations in bold were significant 
in both groups. 
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covariate in the group comparison analyses, differences on simple 
arithmetic word problems and on the coefficient of variation in 
non-symbolic estimation were no longer significant (correspond- 
ing results were also obtained when the influence of the coefficient 
of variation in non-symbolic estimation was removed); and (3) 
contrary to expectations, groups performed equally in the work- 
ing memory tasks. In the following, these results will be discussed 
in further detail. 

COMT AND NON-SYMBOLIC NUMERICAL PROCESSING 

Our most important result is the difference between COMT 
val/val and met+ groups regarding both numerical non-symbolic 
measures (non-symbolic magnitude comparison task and non- 
symbolic magnitude estimation task). This suggests a direct 
link between dopamine regulation and magnitude processing 
mechanisms that is independent from working memory capac- 
ity. Results are better contextualized in terms of the accumula- 
tor model (Meek and Church, 1983). The accumulator model 
assumes that dopamine modulates at least two mechanisms 
involved in non-symbolic number processing: working memory 
and clock speed (for run mode) or counter (for event mode). 
Research has traditionally focused on the role of vall58met 
COMT polymorphism in working memory variability (Mier 
et al., 2010). For the first time, we obtained evidence that COMT- 
related variability, probably related to dopaminergic mecha- 
nisms, may be more directly involved with magnitude processing. 
Variations in dopamine bioavailability alter the clock speed, and 
could affect the precision with which magnitude estimations and 
comparisons are performed. The present results are in line with 
those reported by Wiener et al. (2011), who found that individ- 
uals homozygous for valine allele presented higher coefficients of 
variation in a time comparison task (supra-second) compared to 
individuals with at least one methionine allele. The present results 
also indicate that val homozygous presented poorer acuity in 
the non-symbolic magnitude comparison task, as shown by their 
larger Weber fraction values. Besides that, mathematical learning 
disabilities can be thought as a model to this relationship between 
time and magnitude. The deficit that these children show in non- 
symbolic number tasks has already been well-described (Piazza 
et al, 2010; Costa et al, 2011; Mazzocco et al., 2011; Mejias 
et al., 2011). Furthermore, a recent paper has shown that chil- 
dren with mathematical learning disabilities also present poorer 
temporal estimation when compared to typical peers (Hurks and 
van Loosbroek, 2012). Therefore, both time and quantity modal- 
ities seem to be associated to specific COMT polymorphisms, as 
predicted by the accumulator model (Meek and Church, 1983). 

The present results also suggest that children in the met+ 
group had a more accurate non-symbolic magnitude represen- 
tations. Val/val carriers had higher Weber fraction values which 
means that they are less able to discriminate accurately between 
magnitudes. In a similar way, children with learning disability 
also showed higher Weber fraction values (Piazza et al., 2010; 
Costa et al., 2011; Mazzocco et al., 2011). In contrast, better per- 
formance on non-symbolic numerical tasks has been associated 
with better performance in mathematical tasks in previous stud- 
ies (Halberda et al., 2008). We found that Weber fraction has an 
impact in arithmetic word problems but not in the transcoding 



tasks (Arabic number reading task and Arabic number writing 
task): the more accurate in their non-symbolic skills the subject 
is, the better is the performance on the arithmetic word prob- 
lems. These results are in line with the Barrouillet et al. (2004) that 
claims that transcoding abilities are asemantic, which means they 
do not require access to magnitude representations. However, this 
does not happen with the arithmetic word problems, which com- 
prise only single-digit numbers and probably automatically access 
the magnitudes they represent. Available evidence shows that in 
early stages of development, when they are learning to perform 
single digit operations, children extensively activate prefrontal 
and intraparietal areas. As proficiency is acquired and operations 
are progressively automatically solved as arithmetic facts, the acti- 
vation focus shifts from the prefrontal and intraparietal sulcus 
to the angular gyrus (Rivera et al., 2005). Intraparietal activation 
levels during single-digit arithmetic are also predictive of math 
performance in high school students (Price et al., 2013). 

The similar results obtained when we controlled for the influ- 
ence of both non-symbolic magnitude comparison and esti- 
mation suggest that these two skills interchangeably influence 
performance on arithmetic problems, but not on transcoding 
tasks. This indicates that these two measures tap a common non- 
symbolic (analogic) construct. However, as correlations between 
Weber fraction and the coefficient of variation in estimation were 
significant but of low magnitude (Table 3), other sources of vari- 
ance may influence performance in these tasks. Although both 
tasks are used to access non-symbolic magnitudes, they may also 
recruit other cognitive components (Gilmore et al, 2011). One 
difference could be related to the connection to symbolic repre- 
sentations, more salient in the estimation task. In the context of 
developmental dyscalculia, Rousselle and Noel (2007; Noel and 
Rousselle, 2011) have suggested a distinction between deficits in 
non-symbolic representations and deficits in symbolic to non- 
symbolic access or mapping mechanisms. Performance in these 
two kinds of tasks is, for example, dissociable in children with 
math learning disability (see also Mejias et al., 2011). Some 
authors claim that children with developmental dyscalculia expe- 
rience more difficulties in the symbolic than in the non-symbolic 
tasks (Rousselle and Noel, 2007; Mejias et al., 2011). In the case of 
our results, this could mean that val-val individuals' performance 
in the estimation task improves comparatively less in comparison 
to met+ individuals because it is a more complex task due to its 
symbolic requirements. 

COMT AND WORKING MEMORY 

Although some studies describe differences in working mem- 
ory tasks between groups of COMT polymorphisms (Goldberg 
et al., 2003; Mattay et al., 2003), these results are not consistently 
reported in the literature (Tan et al., 2007; Blanchard et al., 2011; 
Wardle et al, 2013), especially for children (Diamond et al., 2004; 
Mills et al., 2004; Wahlstrom et al., 2007). Distinct arguments may 
help to explain why working memory performance did not differ 
across COMT polymorphic groups. The first one regards single 
nucleotide polymorphism effects and the different approaches the 
authors organize genetic data (see review in Mier et al., 2010). 
Some studies compare all three genotypes (val/val, val/met, and 
met/met) and other consider a genetic dominance as we did in the 



www.frontiersin.org 



August 2013 | Volume 4 | Article 531 [ 7 



Julio-Costa et al. 



COMT and numerical cognition 



present study (val/val vs. met carriers or met/met vs. val carriers). 
Secondly, tasks vary from one study to the other. In our study we 
use digit span and Corsi blocks, but most of the previous research 
use the n-back paradigm that demands more "storage" capacity 
in contrast of the "manipulation" capacity required by our tasks. 
Bruder et al. (2005) proposed that these two kinds of abilities may 
not involve the same circuitry in PFC. As suggested by Wardle 
et al. (2013), some working memory tasks may be more sensitive 
to COMT effects than other ones. Thirdly, the study that found 
differences between groups using the same tasks as in our study 
investigated an adolescent sample (age ranging between 9 and 17) 
(Wahlstrom et al., 2007) and just Vi of our sample was older than 
10 years old. In this way, the results are not easily comparable. 
Moreover, Dumontheil et al. (2011) showed that the difference 
favoring met allele emerges only after 10 years of age, what can 
easily be conciliated with our results. At last, also consistent with 
our results is the finding reported in many studies, which were 
not able to observe behavioral effects in working memory tasks, 
but obtained significant differences in brain activation patterns 
of individuals with COMT polymorphisms (Meyer-Lindenberg 
et al, 2006; Caldii et al, 2007; Tan et al, 2007). 

THE INTERPLAY OF COMT, WORKING MEMORY, AND SYMBOLIC 
NUMERICAL PROCESSING 

In the present study, differences between polymorphisms arose 
in transcoding tasks, which knowingly require working memory 
resources. It is one of the cognitive mechanisms most consis- 
tently associated with Arabic number writing (Barrouillet et al., 
2004; Camos, 2008). According to Barrouillet et al. (2004), ver- 
bal working memory plays a central role in the management of 
the slots to be filled during the transcoding procedure and this 
relationship is even more pronounced in more complex num- 
bers (such as 4 digit-numbers with an internal zero, e.g., 6047). 
Camos (2008) also claims that working memory load is crucially 
involved in the acquisition of more complex transcoding rules. In 
her study, it was found that children with higher working mem- 
ory span presented less transcoding errors, especially the ones 
related to the knowledgment of place-value concept. Even though 
we did not find group differences in our working memory tasks, 
one can assume that underlying differences on this skill between 
COMT polymorphisms might play a role in the variability on 



the transcoding tasks. One possible explanation is that some pre- 
vious studies that found association between working memory 
and transcoding skills used different measures, such as dual tasks 
(Camos, 2008), while our measures of working memory are asso- 
ciated with storage and manipulation of verbal as well as spatial 
information. Zuber et al. (2009) found an association between 
Corsi blocks and number writing performance. That study was 
conducted in German speaking children and thus this association 
is possibly due to the spatial demands related to the inversion rule 
in German language. 

In the arithmetic word problems task, we also found dif- 
ferences between groups. Methionine carriers children had a 
better performance in the arithmetic word problem task. As we 
did not find polymorphic group differences regarding working 
memory tasks, the group effect on arithmetic word problems 
could be explained in terms of differing activation patterns of 
non-symbolic magnitude representations as discussed above. 

In conclusion, these findings suggest that dopaminergic mech- 
anisms may be implicated in numerical magnitude processing. 
Significant differences between children with val/val vs. met+ 
genotypes were also detected on aspects of arithmetic achieve- 
ment that demand more non-symbolic magnitude processing, 
such as single-digit arithmetic word problems. Our results suggest 
that COMT- related dopaminergic mechanisms are implicated not 
only on working memory- related arithmetic performance, as pre- 
viously shown (Tan et al., 2007), but also on non-symbolic mag- 
nitude processing. Further investigation is necessary to clarify the 
association between COMT polymorphisms, working memory, 
and numerical cognition. 
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